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Introduction: A vast, global network of more than 
3200 lobate thrust fault scarps has been revealed in high 
resolution Lunar Reconnaissance Orbiter Camera 
(LROC) images [1-6]. The fault scarps very young, 
<50 Ma, based on their small scale and crisp appear- 
ance, crosscutting relations with small-diameter impact 
craters, and rates of infilling of associated small, shal- 
low graben and may be actively forming today [1-3]. 
The population of young thrust fault scarps provides a 
window into the recent stress state of the Moon and 
offers insight into the origin of global lunar stresses. 
The distribution of orienations of the fault scarps is non- 
random, inconsistent with isotropic stresses from late- 
stage global contraction as the sole source of stress 
Modeling shows that tidal stresses contribute signifi- 
cantly to the current stress state of the lunar crust [1]. 
Tidal stresses (orbital recession and diurnal tides) su- 
perimposed on stresses from global contraction result in 
non-isotropic compressional stress and thrust faults con- 
sistent with lobate scarp orientations. Stresses due to 
orbital recession do not change with orbital position, 
thus it is with the addition of diurnal stresses that peak 
stresses are reached. At apogee, diurnal and recession 
stresses are most compressive near the tidal axis, while 
at perigee they are most compressive 90° away from the 
tidal axis [7]. Coseismic slip events on currently active 
thrust faults are expected to be triggered when peak 
stresses are reached. Analysis of the timing of the 28 the 
shallow moonquakes recorded by the Apollo seismic 
network shows that 19 indeed occur when the Moon is 
closer to apogee, while only 9 shallow events occur 
when the Moon is closer to perigee [7]. Here we show 
the results of relocating the shallow moonquake using 
an algorithm designed for sparse networks [12] to better 
contrain their epicentral locations in order to compare 
them with stress models. The model for the current 
stress state of the Moon is refined by investigating the 
contribution of polar wander. 
Current Stress State and Polar Wander: Radial con- 
traction from interior cooling is the dominant source of 
stress, contributing >2 but <10 MPa based on the cur- 
rently mapped population of lobate scarps [1, 2]. Super- 
imposed on compressional stresses from contraction 6- 
are two components of tidal stress, orbital recession 
stress o, and diurnal stress o;. Tidal stresses are domi- 
nated by o; that may reach 20 to 40 kPa [1]. An addi- 
tional component of stress that may significantly con- 


tribute to the current lunar stress state is polar wander. 
True polar wander has been attributed to a change in the 
Moon’s moments of inertia due to a low-density thermal 
anomaly beneath Procellarum [8]. The change in the 
location of the poles is consistent with the observed 
remnant polar hydrogen deposits [8]. Modeling shows 
that stresses from ~3° of polar wander oy, over the last 1 
billion years results in stresses with orientations and 
magnitudes similar to o; (Fig. 2). The net non-isotropic 
compressional stresses Go, from 6;+Oyw+o- result in thrust 
faults with preferred orientations that are in general 
agreement with orientations of the mapped faults (Fig. 
3). Although the contribution of diurnal tidal stresses to 
On is small (<5 kPa), the addition of o; results in peak 
stress when the Moon is near apogee or perigee, con- 
sistent with the occurrence of most shallow seismic 
events [7]. 

Shallow Moonquakes: Four seismometers were placed 
on the Moon at the Apollo 12, 14, 15 and 16 landing 
sites. During the operation of the seismometers (1969 to 
1977), 28 shallow moonquakes were recorded [9, 10] 
(Fig. 1). Shallow moonquakes occur at depths <200 km 
and have been interpreted as tectonic in origin [10]. 
Analysis of all 28 shallow moonquakes indicates Rich- 
ter-equivalent magnitudes in the range from 1.6 to 4.2 
[10]. Estimates of stress drops range from a few MPa to 
over 100 MPa for the 3 largest events [11]. 

The current best estimates for most shallow moon- 
quake locations are likely only accurate to several de- 
grees, possibly more for small and/or distant events. 
This makes a comparison with stress models and at- 
tempts to match the exact locations of known tectonic 
features with specific shallow moonquakes problematic. 
Their depths are similarly ill-constrained, with some 
estimated to occur at the surface, and others up to 200 
km deep with uncertainty up to an additional 200 km. 

The standard method for locating a seismic event uses 
a known velocity model and the observed arrival times 
of the direct P and S waves. Large uncertainties in arri- 
val times directly translate into a large uncertainty in 
event location. We have instead applied a location algo- 
rithm specifically designed for sparse networks [12] to 
determine whether their location accuracy can be im- 
proved. Rather than solving for a best-fit location, this 
approach divides the solution set into falsified and non- 
falsified candidate locations using an adaptive grid 
search, and accounts for arrival time uncertainty using 


windows around the true arrival time. Preliminary re- 
sults suggest that as many of half of the published shal- 
low moonquake locations may fall outside the non- 
falsified regions. 
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Figure 1. Map of lobate scarps (red), epicentral locations of shallow moonquakes (blue dots), and locations of Apollo Seismic Network seismo- 
meters (black diamonds). Moonquakes are scaled by estimated Richter magnitude [10]. Mare basalt units are shown in tan. 
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Figure 2. Orientations of principal stresses due to orbital recession and polar wander. Compressional stress components are show by red lines and 
reach a maximum of ?? kPa. Extensional stress components are shown by blue lines. 
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Figure 3. Agreement between observed faults and expected local thrust fault orientations due to different stress models. Fault segments at 0 
on plot are orthogonal to local direction of most compressive stress, and at 90 they are parallel to most compressive stress. 


